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COMPUTER PROGRAM TO SIMULATE RAMAN SCATTERING 


By 

Barbara Zilles 1 and Roscoe Carter 2 


SUMMARY 

A computer program is being developed to simulate the vibration 
rotation and pure rotational spectrum of a combustion system con- 
sisting of various diatomic molecules and CO 2 as a function 
of temperature and number density. The main program for this 
purpose has been developed along with the necessary subroutines to 
generate two kinds of spectra — a pure rotational spectrum for any 
mixture of diatomic and linear triatomic molecules, and a vibra- 
tional spectrum for diatomic molecules. The program is designed 
to accept independent rotational and vibrational temperatures for 
each molecule, as well as number densities. 

In the following section, the theory upon which the program is 
based is described in detail. In the third section block diagrams 
of the program are given. The different variables and functions are 
identified and their relation to the theory and/or logical flow of 
the computing format is explained. In the "results" section example 
of computed spectra are given along with a brief description of the 
input parameters used to produce these spectra. 


1 Research Staff, Old Dominion University Research Foundation, 
Norfolk, Virginia 23508. 

2 Assistant Professor of Chemical Sciences, Old Dominion University 
Norfolk, Virginia 23508. 



THEORY 


The power signal collected for a Raman transition from rotation- 
vibration state Jv to J 1 v ' is given by 


R 


J' v' 

Jv 


n Y jv Vlfte 



( 1 ) 


where n 


Y 


Jv 


V 

I 

e 

a 


= number density 

=s probability of finding a molecule in rotation- 
vibration state Jv according to the Maxwell- 
Boltzmann distribution law 

= irradiated volume 

= average intensity of irradiation in that volume 
= efficiency of collecting optics 
— scattering angle in steradians 



J' v' 


Jv 


differential Raman cross section for the transition. 


Optical Factor 

The problem of defining a volume in the sample which interacts 
with the laser light and the average intensity within that volume 
has been considered by Barrett and Adams (ref. 1). They have 
shown that for a laser beam focused to a diffraction limited point 
by a lens of focal length equal to that of the collecting optics 
a cylindrical volume may be defined 

V = 2bA f , (2 

where A^ is the area of the laser focal point and b is the 
distance in both directions along the beam axis at which the 
irradiance has dropped to 1/e 2 its intial value. A f is 
related to b by 


2 



(3) 


b = 2A f /A 


where A is the wavelength of the laser line. The average intensity 
in this volume is found to be 


I = 



tan~ 1 (2) I Q 
8 


(4) 


where A is the area and I the irradiance of the unfocused laser 
u o 

beam. Substituting equation (3) into equation (2) and combining 
with equation (4) gives 


VI = 


A A £ 
u f 

2 A 


jtan 1 (2)J 1^ 


Now, A f is related to A^ by 

A r = 4A 2 f 2 /A 
f u 

where f is the focal length of the focusing lens. Hence we have 

VI = 2Xf 2 |tan~ 1 (2)1 I o • < 5 > 

The irradiance I is given by 

I - P/A (6) 

o u 

where P is the laser power. 

Frequency 

The frequency associated with the transition giving rise to the 
signal in equation (1) may be separated into rotational and vibra- 
tional terms. 
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( 7 ) 




J'v 

Jv 


= to 


V 
3 

V 


+ 03 


= G ( V ' )l - G(v) + [f v ,(J') - F v (J) 


where G and F are the vibrational and rotational energies, 

v 

respectively, of the states in question. The rotational energy 
depends on the vibrational as well as the rotational quantum 
numbers : 


F (J) = B J (J + 1) 
v v 


( 8 ) 


If B q is the rotational constant of the ground vibrational state, 


b = b -y a. 

v 0^1 


. V. 
1 


(9) 


where a. is the vibration-rotation interaction constant and v. 

i l 

is the vibrational quantum number for the ith normal mode of 
vibration. 

For the pure rotational spectrum, of course, v f - v and 

XT I .... 

to =0. For vibrational spectra the transitions will involve 
a change of quantum number for only one normal mode (that is, 
combination bands will not be considered) . Thus, dropping the 
normal mode subscript for Av = v' - v and a, the rotational 
part of the transition frequencies for the different branches of 
the band are as given in table 1. [B^ in the table, however, 

still denotes a sum over normal modes as in equation (9)]. 

The vibrational energy relative to that of the ground state 
is given by 




o 

x . . 

ID J 


v . 
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-E 

i 


£ 

j<i 



v . v . 
i D 


( 10 ) 
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j. (a) 

Branch J' U) T 

J 

0 J-2 -4B v (J - - a (J - 1) (J - 2 ) Av 

P J-l _2B V J “ aJ(J - 1) Av 

Q J -a J ( J + 1) Av 


R j J+l j 2B v (J +1) - a(J + 1) (J + 2 ) Av 
4B v (J + -|) - a(J + 2) (J + 3 ) Av 


S 


J+2 



S £ 

(b) 



JJ 1 



3 (J 2 

-■a 2 ) [(J - D 2 

- £ 2 ] 

2 J ( J 

- 1) (2 J 

- 1) (2J 

+ 1 ) 


3£ 2 {J 2 

- I 2 ) 


J(J + 

i) (j - 

1) { 2 J - 

u 

[j(j + i) - 

3£ 2 ) ] 2 


J(J + 

1) ( 2 J - 

1) (2 J + 

"Tf 

3it 2 

[ (J + 1) 

2 - l 2 ] 


J(J + 

i) (j + 

2) (2J + 

IT 


3 [ (J + 1) 2 - £ 2 ] [ (J + 2) 2 - l 2 ] 
2 (J + 1) (J + 2) ( 2 J + 1) ( 2 J + 3) 
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a). is the harmonic frequency of the ith vibrational mode 

x 0 . 1 

ID 

and i. For the band involving transitions of the kth normal 


the anharmonicity constant for the interaction of modes 


where 
and 

iD 

1 and j 

mode, therefore, the vibrational part of the transition frequency 
is given by 


03 


(k) = 03 %v - x° (2v k + Av + l)Av - £ x° k (v ± + ~) Av 

i^k 


where Av denotes the quantum number change, that is, Av = v k - v k « 
For the fundamental band of a diatomic molecule, there is only one 
normal mode and Av = 1. Hence the above equation becomes 


03 V = 03° - 2x° (v + 1) 

v 


(ID 


Boltzmann Factor 

The population of the Jv rotation-vibration state is given 


by 


Y JV = g J (2J + X) e 


-F(J)/kT- -G (v)/kT 


v 


/Z Z° 
r v 


( 12 ) 


where gj is the nuclear degeneracy of the rotational state, T^ 
is the rotational temperature, and T^ is the vibrational tempera- 
ture. F (J) and G° (v) are the rotational energy and the vibra- 
tional energy relative to the zero point energy, respectively. 

Z r is the rotational partition function and Z° is the vibrational 
partition function relative to the zero point population. The 
rotational partition function can be approximated by an integral 
(see for example ref. 2, Appendix C) : 


~ -(B - a . v . ) J (J+l)/kT 

z =£ gj <2J + 1) TTe ° 11 

i 


** \ kT r (g + + g )/E ( B 0 " a i v i) 

i 


(13) 


6 



where g + and g~ are the nuclear degeneracies of the even and 
odd J levels. Thus the partition function product is equal to 


Z Z° = ~ kT (g + + g ) Z(T ) 
r v 2 r ^ ^ v 


-G°(v)/kT 
g e v 

-E ? J— ^ 

v B -V a . v . 

o ^ 11 


( 14 ) 


where G°(v) is given by equation (10) and g^ is the vibrational 
degeneracy — which equals 1 except for states (occurring in C0 2 ) where 
a degenerate normal mode, k, is excited, g^ = v^ + 1. 


Cross Section 

The differential scattering cross section in equation (1) is 
given by (ref. 3) : 



2 4 tt 4 oo 4 |<J* v* | a | Jv>| 2 

Ja. 


(15) 


to is the scattered Raman frequency, 
R 


03 


R 



03 


J* V* 

Jv 


(16) 


J ^ v ^ 

where A is the exciting wavelength and is the frequency 

of the transition given by equation (7). 

<J , v'|a|Jv> is the matrix element of the polarizability tensor 
for the initial and final vibration- rotation states in the labora- 
tory-fixed system. In terms of the molecule-fixed coordinate system 
the diagonal tensor element (scattered light polarized parallel to 
the incident electric vector) is given by 


a 


ZZ 



cos 


(if Z) 


CL. . COS 
1] 


(j f Z) 


(17) 
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For linear and symmetric top molecules, if ou ^ is broken down 
into a sum of spherical and traceless tensors (ref. 3) 

|<J'v' |a zz |Jv>| 2 = I<v'|a|v>| 2 6 JJt 

(18) 

+ 43 l<v’ I 3 I v>| 2 Sjj, 


where <5__, is the Kronecker delta and S , is the rotational 
line strength. The parameters, a and 3, referred to as the 
polarizability and anisotropy, respectively, are two invariants 
of the polarizability tensor a 



They are functions of the normal vibration coordinates Q and may 
be developed in a power series 


a (Q) 



Q + . . . 


3(Q) 


o 



Q + 


/ 


Thus 


|<^'|a|v>| 2 = a 2 6 vv , + a’ 2 |<v + l|Q|v>| 2 <5 v , , v+1 (20) 

where the zero subscript has been dropped. For a diatomic molecule 
the normal coordinate matrix element is given by 
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( 21 ) 


| <v + 1 | Q | V> |2 = Nh(v i) , 

8tt 2 cco y 

where N = Avogadro ' s number 
h = Planck's constant 
c = speed of light 

w 0 = the harmonic frequency as in equations (10) and (11) 
y = the reduced mass. 

Similarly, the anisotropy for any molecule may be expanded 

|<v'|3|v>| 2 = 3 2 6 vv , + 3 1 2 |<v + l|Q|v>| 2 6 v ., v+1 ( 22 ) 

The first term in equation (22) when substituted into the second 
term in equation (18) gives the pure rotational spectrum for parallel 
polarization. The first term in equation (20) gives the contribution 
to the Rayleigh scattering. 

The second terms in equations (20) and (22) determine the 
vibrational spectrum. For perpendicular polarization the off- 
diagonal polarizability tensor element a Y 2 i n the laboratory 
system is considered using a transformation analogous to equation 
(17). For linear and symmetric top molecules, 


|<J'v' |a YZ | Jv>| 2 = A, |<y'|6|v>| 2 Sjj, . (23) 

Thus for vibrational spectra we can define trace and anisotropic 
scattering cross sections: 


da\ 

dft I 
/trace 



— d a \ 

3 dtt 


<v' | a | v>| 2 6 JJt 


(24) 


da\ 

dft j 


anis 


45 da \ 

3 (I 


|<v'|3|v>| 2 




(25) 
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where the parallel and perpendicular cross sections are determined 
by whether the component of a in equation (15) is ZZ or YZ, 
respectively. 


For the diatomic molecules in this study, the parallel scattering 
contains significant contributions from the anisotropic as well as 


the trace cross sections. For polyatomic molecules 



depends on the normal coordinate being considered. It is always 
greatest for the totally symmetric vibration. For spherical top 
molecules, 3^ = 0 for the strong totally symmetric vibration 
and thus only the trace scattering need be calculated. 


For asymmetric top molecules the anisotropic scattering does 
not assume such a simple form as equation (25) . However, for H 2 0, 
the asymmetric top considered in this study, the anisotropic scatter- 
ing is weak compared to the trace scattering and does not contribute 
to the appearance of the spectrum [the depolarization ratio = 1 ^/ 1 ^ 
is 0.02 (ref. 4)]. The experimental trace scattering may be deter- 
mined quantitatively since the relation in equation (24) still 
applies. Hence for H 2 0, only the trace scattering will be calculated. 


The rotational line strength factors S , are determined by 

O J 

the matrix elements of the direction cosines as seen from equations 
(17) and (18). They determine the rotational selection rules. For 
pure rotational spectra the selection rules are 


AJ = J* - J = ±2 ' (26) 

for the diatomic molecules. Since the antistokes scattering (AJ = -2) 
is symmetrical to the stokes scattering (AJ = +2) , except for the 
factor in equation (15) and the Boltzmann factor, it gives no new 
information and is not calculated. 

For the triatomic linear C0 2 molecule, those states where the 
doubly degenerate bending mode, . v 3 , is excited, can have non-zero 
vibrational angular momentum quantized along the figure axis. This 
angular momentum vector contributes to the direction cosine matrix 
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elements in the same way as the quantized component of the rotational 
angular momentum for a symmetric top. Thus for states with a non- zero 
vibrational angular momentum quantum number Z the line strengths 
depend on Z as well as J and J' and we have the selection 
rules 


AJ = ±1 for Z ^ 0 (27) 

in addition to equation (26) . 

For C0 2 states with v 3 0, there are v 3 + 1 states with 
the same energy. When v 3 is odd all the states have Z ^ 0 and 
there are two states each with | it | = 1, 3, 5, # .«v 3 . When v 3 is 
even there is one state with Z = 0 and two states each with 
| % [ = 2, 4, 6 ,...V 3 . For the pairs of states with the same Z value, 
a given J level of only one of the states is occupied for C0 2 , 
since one state is symmetric with respect to interchange of the 
identical oxygen nuclei and one state is antisymmetric. Since 
oxygen has zero nuclear spin, only the even J levels are occupied 
for symmetric states and only the odd J levels for antisymmetric 
states. States with Z = 0 are symmetric. 

In addition to different nuclear degeneracies, states with 

different Z have different line strength factors so that S' 

J J 

in equation (18) is given by 

s jj- -£ 4 s jj- • < 28 > 

z 

The S , are given in table 1 for AJ = 0, ±1, ±2; they reduce to 
the line strength factors for diatomic molecules when Z = 0. For 
rotational spectra J* = J + 1, J+2; for vibrational spectra all 
five branches are calculated. 

Units and Summary 

J ^ V ^ j 

In equation (1) R is in watts, since VI is in watt-cm, 

J V 

n is in cm "* 3 and the cross section is in cm 2 . We may obtain inten- 
sity in photons/sec by dividing by 
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1 photon at frequency oj_ = hcto^ watt-sec 

K K 

where h is in joule-sec, c is in cm/sec, and is given by 

equation (16) . 

In summary then, for rotational spectra the signal in counts/sec 
is 


o 


R 


J' v* Cn 


Jv 


Z “R g j (2J + 1>e 


-F(J)kT -G (v)/kT 


V 




(29) 


where Z = Z Z /k and 
r v 


c = 2 4 tt 4 4 z 3 2 

4 5 khc 


(30) 


For vibrational spectra of diatomic molecules the 0 and S branch 
signals are 

C 1 n 3 A (v + 1) 

— i - 

Z R o 
a) y 

(31) 

where A = Nh/87r 2 c and C* is the same as C except that 3 is 
replaced by 3'« 

The signal for the vibrational Q branch is the same as equation 
(31) except that s jj i is replaced by 

11.25 a ,2 /3' 2 + S TT . 

J J 


o 


gj (2J 


+ 1) e 


-F(J)/kT -G (v)kT 



R 


J'v 

Jv 


PROGRAMMING 

The computer program consists of a main program to take care 
of the major input/output operations, calculate general constants 
and manage the general flow of operations; a subroutine, PART, to 
calculate the partition functions; a subroutine, LINES, to generate 
the corresponding frequency and intensity arrays; a subroutine, CNVLUT, 
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to convolute with the line shape function; a subroutine, SLIT, to 
convolute with the slit function; and a subroutine, SPKPLT, to plot 
the resulting spectrum. 


The program is set up so that for a given run only a given type 
of spectra may be computed. The alternatives are; 


1. 

a set of 
diatomic 

pure rotational spectra for any combination of 
or linear triatomic molecules. 

2. 

a 

set 

of 

vibrational 

spectra 

for 

diatomic molecules. 

3. 

a 

set 

of 

vibrational 

spectra 

for 

CO 2 t 

4. 

a 

set 

of 

Vibrational 

spectra 

for 

h 2 o. 

5. 

a 

set 

of 

vibrational 

spectra 

for 

ch 4 . 


Thus a different LINES subroutine is placed in the card deck according 
to the desired option. At present the LINES subroutine has been 
completed only for the first two options. 


Main Program 


A flow chart of the main program is shown in figure 1. The 
input parameters and the constants calculated in the main program 
are as follows: 


Instrument parameters. A, exciting wavelength in A; P, laser 
power in watts; A^, area of unfocused laser beam in cm 2 ; ft, 
collecting angle in steradians; f, focal length of focusing lens 
in cm; e, total efficiency. 


Optical factor. F 
equation (5) . 


2 A 


tan 1 



x 10~ 8 f 2 ftPe/A Jl . 


See 


Spectrum parameters. DX, plotting increment in cm" 1 ; Au), 
halfwidth of the line shape function in cm" 1 ; RLIM, limit of 
intensity in counts/sec, below which R for a given line is set 
equal to zero; SLT, spectral halfwidth of the slit in cm" 1 ; 
ALNGTH, length of the plotted spectrum in inches; NSTRT , fre- 
quency in cm" 1 at which spectrum begins; IEND, frequency in cm" 1 
at which spectrum ends. 
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NPT is the number of plotting points. NPT = (NSTRT - IEND) / 

DX + 1. 

Molecular constants. B , ground state rotational constant in 
cm" 1 ; 3 or 3 '/ anisotropy in cm 3 for pure rotational spectra or 

anisotropy derivative in cm 2 for vibrational spectra; ATB, (q’/P 1 ) 2 / 
zero for a pure rotational spectrum; ]i, reduced mass--zero for a 
pure rotational spectrum; NQ, number of normal coordinates — must 
be 1 or 3 ; g and g , nuclear degeneracy of positive and nega- 
tive J levels, respectively; NJ, maximum J quantum number; 
SPECIE, name of molecule. 

= rotational temperature; = vibrational temperature; 
n = number density in molecules/cm 3 . 

Normal mode, Q, parameters. U)°(Q^), harmonic vibration 
frequency in cm- 1 ; x°(Q^), anharmonicity constant in cm"* 1 ; a(Q^), 

rotation-vibration interaction constant in cm"* 1 ; NS(Q^) = +1 or -1, 
symmetry of normal coordinate with respect to reflection through a 
plane containing the figure axis. If the symmetry of the vibra- 
tional wave function is negative, the nuclear degeneracies of the 
odd and even J levels are interchanged. This symmetry is given by 

n =TTnS(Q.) 1 . (32) 

S . 1 

1 


NVM(i), maximum v quantum number for the ith vibrational mode, 
co 1 is a factor used in the PART and LINES subroutines, 

l 


i o o 

(oi = 00. - X. . 

1111 


(33) 


D 


,44 4 P 2 nF 

45 khc Z 


(34) 


See equations (29) and (30) . 

The index L in the main program serves to manage the generation 
and convolution of the different branches of a vibrational spectrum. 

In the rotational LINES subroutine the R branch, if allowed, is 
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calculated simultaneously with the S branch; hence it sets L to 
zero so that at the L flag in the main program the flow will be 
recycled to begin generation of lines for another molecule if 
desired or to plot the spectrum. 

Similarly, if only the trace scattering of the vibrational 
spectrum is to be generated for H 2 0 (because the anisotropic 
scattering is impracticable) or for CH 4 (because the anisotropy 
is zero) , L will be set to zero in the LINES subroutine and only 
the Q branch will be generated. For vibrational spectra in which 
the anisotropic as well as trace scattering is to be computed, 

L = 1, 2, 3 correspond to the 0, Q, and S branches, respectively. 

The index MOL is introduced primarily for the case of the 
pure rotational spectrum where the scattering from several molecules 
appears in the same spectral region. It serves to index the specie 
concentrations and temperatures for labeling of the plotted spectrum. 

The addition of the background, BKGRD, to the array of spectral 
ordinates, Y, has not yet been implemented as this data is not 
available. It is presently assumed that the data will be decon- 
voluted and hence the background is added between the line shape 
and slit function convolutions. 

PART Subroutine 

The subroutine PART calculates Z, the rotational times the 
vibrational partition function divided by the Boltzmann constant k. 

A flow chart of this subroutine is given in figure 2. The variables 
calculated are as follows: 


num = numerator of Z (Ty) without g^ [see equation (14)] 


num 


= TT 



x?. v.)v./kT 
1111 v 


(35) 
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Figure 2 . Flow chart for subroutine PART 














0 for i ^ j , 


as in equation (10) where it is assumed that x?. = 
and w l is given by equation (33) . 

den = denominator of Z (T^) den = B q ou 

i 


Z(T V ) =E 

V 


num(v) (v 3 + 1) 
den (v) 


Z = 


Z Z° 
r v 


+ 


= T 


(g + g ) 


Z(T v ) 


The sums and products in this function are performed numerically. 
Theoretically the sum over vibrational levels is infinite. In prac- 
tice NVM(I), the maximum v quantum number, is input so that 
the population of this level is less than 10 ~ 6 that of the ground 
state. NZTV(I) is then set equal to NVM(I) +3 and higher con- 
tributions to the sum should be negligible. If there is only one 
normal coordinate (diatomic case) NZTV(2) = NZTV(3) = 0. 

Rotational LINES Subroutine 

The rotational LINES subroutine calculates the frequencies and 
intensities of the spectral lines for the R (C0 2 ) and S branches 
according to equation (29) and table 1. A flow chart of this sub- 
routine appears in figure 3. Its rather complicated appearance is 
due to the accomodation of linear triatomic molecules with their 
three degrees of vibrational freedom and non-zero vibrational 
angular momentum. The loops over quantum numbers for the different 
vibrational modes are nested in order of decreasing frequency 

> a) 2 > 0 ) 3 . The loop over J is nested within the loops over 
v, and the loops for summing over the degenerate vibrational 
substates (£ quantum number) as in equation (28) are nested 
within the J loop. The parameters and functions calculated 
are as follows. 
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NVM (2) •--- 0 


| NVM ( 3) = 


calc to, a 
for S branch 


cnl.c R 
for S branch 


NVTOT = 0 


DOX^V! - 0 , NVM ( 1 ) 

1 ~Z 

DO\^7v 2 = 0 , NVM ( 2 ) 


NVTOT = NVTOT + 1 


RS (NVTOT, S) = 0 


DO^Vg = 0 , NVM ( 3 ) 


DO I = 1 , NQ 


calc to,S 
for R branch 



calc R 
for R branch 


NVTOT = NVTOT + 1 


T = v 3 .EQ.O 


DO\o/j = 0,NJ 


R (NVTOT , J ) = 0 ,NJ 


I = n (-1)' 



i = £ Q ,v 3 , 2 




g T - g + g 


Figure 3. Flow chart for rotational LINES subroutine 
























NVTOT is the index of vibrational states for all the different 

, ^ J+2,v 

v lf v 2 , v 3 and also for the R and S branches. Thus if o) Jv 
is given by tx> (NVTOT, J) , o)^ 1,V is given by w (NVTOT+1, J) . Thus 

NVTOT must be decremented for each cycle over J if v 3 ^ 0 and 
must be re-incremented before initiating a new cycle over v 3 . 

B is given by equation (9) and n is given by equation (32) . 
v s 

E is the exponent in the vibrational part of the Boltzmann 

factor [see equation (35)], 


E 


= L 

j<i 


(UJ* - 

3 


x 


33 


V . ) V . 

3 3 


If the above sum exceeds 13 for a given v 3 , then any states with 
the same v 1 ,v 2 values will have even higher populations for higher 
v 3 , so control is transferred out of the v 3 loop; similarly for 
v 2 . For Vj control is transferred to the end of its loop. 

I = l for Vq odd; £ =0 for v 3 even as in the discussion 

o ° o 

following equation (27) . 

T is a logical variable, tested for calculating R branches 
and processing NVTOT. 

+ — , 

g is the nuclear degeneracy set equal to g or g according 
to whether the vibration-rotation product wavefunction is symmetric 
or antisymmetric. After the contribution to R for the £ = 0 state 
for a given v 3 has been calculated, gj is set equal to g + + g 
since £ ^ o states occur in pairs as discussed following equation (27) . 

BF is the Boltzmann factor without g T and the (2J+1) factor 

J £ 

which cancels with that in the denominator of S jjt for the R and 
S branches (see table 1) . 

(J.LT.£) is a logical variable which if true transfers control 
to the end of the J loop, since J is the total angular momentum 
quantum number and states with J < £ cannot exist. 

to and S , are as given in table 1 with Av = 0, except that 
the ( 2 J+l) factor in the denominator of S* » has been cancelled out. 
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R is given by equation (29) using the factors which have 
been calculated: 

R (NVTOT, J) =£ D co(NVTOT,J)^ g^(BF)Sj Jf 

where D was calculated in the main program according to equation (34) . 

Vibrational LINES Subroutine 

The vibrational LINES subroutine calculates the frequency and 
intensity of the spectral lines for the 0, Q, and S branches 
of a diatomic molecule according to equation (31) . The flow chart 
is given in figure 4. It follows the same general form as the 
rotational routine without the sums and products over multiple 
vibrational coordinates and allowing for the difference in formulas. 

The variables calculated are as follows. 

DD is the constant factor D in equation (34) weighted by 
the normal coordinate matrix element, as given by equation (21) , 
without the (v+1) factor, 

co_ is the vibrational part of the transition frequency for 
F 

v = 0, u)_ = w 1 - x° according to equations. (11) and (33). 

F 

IEX is a factor used in calculating c 2 (see below) , 

IEX = -1 for L = 1; IEX - 1 for L = 3. 

Ci is the coefficient for the J independent terms in the 
frequency formulas for the 0 and S branches, c x = - a. 

co is the vibrational part of the transition frequency, 

~ o 

co = (o_ - 2x v. 
v F 

co includes the J independent terms from the rotational 
c 

part of the transition frequency, co = w + 20 ^ . 

c 2 is the coefficient of the frequency terms linear in J, 
c 2 = 4cj (IEX) - a . 

H is the DD factor weighted by (v+1) — see above. 
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(J.LT.2) is a logical variable tested to eliminate negative 
J" values for the 0 branch (J' = J + 2) . 

g and BF are the same as in the rotational subroutine. 

J 

T * 

oo and S , are as given in table 1, with i = 0 and Av = 1, 
J J J 

00 = 00 + CoJ - aJ 2 for the 0 and S branches. And oo = oo - 

c ^ v 

aJ(J+l) for the Q branch. The (2J+1) factor in denominator of 

S | for the O and S branches has been cancelled out with 
J J 

that in BF as in the rotational subroutine. 

R is given by equation (31) using the factors which have been 
calculated. Thus for the 0 and S branches. 


/ 1 3 

R (NVTOT , J) = H ( j - w(NVTOT,J)^ g a (BF)S JJf . 

And for the Q branch, 

R {NVTOT, J) = H (2J+1) - w (NVTOT, J)^ 3 gj (BF) (11.25 ATB + Sjj) . 

where ATB = (a'/3') 2 « 


CNVLUT Subroutine 

The CNVLUT subroutine calculates the ordinates Y (N) of the, 
spectrum by convoluting over the generated lines with a Lorentzian 
line shape function. The flow chart is given in figure 5. 

(MOL.EQ.l) and (L.GE.2) are flags which insure that the 
cumulated data are saved from a previous molecule (pure rotational 
spectra) or from previously convoluted branches (vibrational spectra) . 
Otherwise the array of ordinates, Y (N) , is initialized to zero. 

YMAX saves the maximum ordinate value for scaling the spectrum 
before plotting and labeling it. 

X is the abscissa in cm -1 . 

(J.LT.2) has the same function as in the vibrational LINES 
subroutine. 
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Figure 5. Flow chart for subroutine CNVLUT. 
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The Y (N) values for each point are calculated according to 


*<n> = 

v J 


R(v, J) (Aw) 2 

(Aw) 2 + [w (v, J) - X] 2 


The C and RLIM flags are to eliminate unnecessary calcu- 
lations in this time-consuming routine, C = X - w(v,J). If 
| C | > 45 the contribution of R(v,J) to Y (N) is neglected. 

If w(v,J) < X (C > 0) , the program skips over J values until 
w is brought within the range of 45 cm” 1 from X, since w 
increases with increasing J values. If w(v,J) > X (C < 0), 
then once |c| > 45 cm -1 no higher J values will bring it into 
this range so control is transferred outside the J loop to the 
next v — w(v,J) decreases with increasing v. Actually for 
the O and Q branches w decreases with increasing J (see 

table 1). The (L.EQ.l) or (L.EQ.2) flag, therefore, sets 
C = -C so that the C flags will have the same effect for 
the O and Q branches as for the R and S branches. 


The (X - w„) flag allows for truncating the Q branch of 
a vibrational spectrum so that the structure of the rotational wings 
can be observed. Thus an ordinate is not considered maximum unless 
it is 10 or more cm” 1 from the v = 0 line (w„) of the Q branch. 

r 

In the rotational LINES subroutine w p is set equal to -20 so 
that (X - w„) will never be less than 10. 

r 

SLIT and SPKPLT Subroutines 

The SLIT subroutine is set up analogous to CNVLUT except that 
the convolution function is triangular rather than Lorentzian and 
the sum is over the data points generated by CNVLUT rather than the 
quantum numbers. It also generates and stores an array of abscissa 
values, x (N) , corresponding to Y (N) . At present there seems to 
be something wrong with this subroutine as the program has never 
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completed a run since it was added. The DEC-10 computer throws the 
job off the system without leaving a LOG: FILE or any record of the 
job having been logged on. 

The SPKPLT subroutine generates a PLOT: FILE from the ordinate 
and abscissa arrays. It consists primarily of subroutine calls 
which are specific to the DEC-10 system. This subroutine is 
currently being rewritten for use on the CDC-6000 computer and 
is not given here. 


RESULTS 


Examples of computed spectra appear in figures 6 through 14. 
For all these spectra the same instrument parameters were used. 
They are: 


o 


exciting line = 4880 A 

laser power = 0.5 watts 

area of unfocused 

laser beam = 0.0128 cm 2 

collecting angle = 0.05 steradians 

focal length = 50 cm 

efficiency - .10 


The convolution over the triangular slit function was not per- 
formed for these spectra. Only a single convolution with the 
Lorentzian line shape was performed, but the halfwidth used was 
typical of (the lower limit of) spectral slitwidths, that is, 

0.5 cm” 1 — whereas natural line widths would be on the order of 
0.01 to 0.1 cm”” 1 . A typical plotting increment was .2 cm”” 1 . 

RLIM was chosen as 1 x 10~ 6 for the vibrational spectra and 
1 x 10~ 3 for the pure rotational spectra. 

The molecular parameters used in the computed spectra are 
given in table 2 and the normal coordinate parameters in table 3. 
The B , a, u)?, and x? . values were derived from analysis of 
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Figure 7. Rotational Raman spectrum of 78% N 2 , 15% 0 2 , 7% C0 2 at 600 
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Rotational Raman spectrum 
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N 2 , 15% 0 2 , 7% C0 2 
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Figure 13. Vibrational Raman spectrum of 0 at 600 °K. 
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Table 2 


B cm" 1 
o 

3 , 10~ 24 cm 3b 

i 

,10- 16 cm 2b 

ATB b 

+ 

g 

2.00065 a 

r ' - ■ ■ ■ 1 

.950 

2.16 

.546 

6 

1.43776 a 

00 
i — ! 

• 

i— 1 

2.34 

.357 

0 

.39021° 

2.44 



1 


Reference 5. 
Reference 6. 
Reference 7. 


Table 3 


0 

03 . 

1 

o 

x i± 

a . 
i 

1 

NS . 

l 

2359. 61 a 

14.456 3 

,0187 a 

1 

1580. 36 a 

12.073 a 

. 0157 9 a 

1 

2349. 16 b 


.00309 d 

-1 

1388. 17 b 


.00121 d 

1 

667. 4 b 


rt 

-.0072 

1 


Reference 5, 


Reference 8. 


The numbering of normal coordinates is that used in 
the program — 0 ) x > o> 2 > 0 ) 3 — and does not correspond 
to standard numbering. 

. L- 

Reference 7. 

























infrared spectra (refs, 5, 7) . C0 2 anharmonicities were not used 

for the spectra given here; however, these data are available 
and should be used. 

3 was derived from measurements of the depolarization ratio 
of the Rayleigh band (including rotational wings) combined with a 
calculated from refractive index measurements (ref. 6) . ATB was 
derived from the depolarization ratio of the Raman vibration-rotation 
band, and 3' was derived from the latter combined with the inten- 
sity ratio of the Raman to Rayleigh bands (ref. 6) . 

NVM and NJ, the maximum v and J quantum numbers, are 
given in table 4. They have been chosen to correspond to the level 
whose population is less than 10” 6 that of the most populated level, 
for the given temperature. 


Table 4. 



NVM 

NJ 

298° 

600° 

1000° 

298° 

600° 

1000° 

n 2 

2 

3 

5 

40 

58 

76 

o 2 

2 

4 

7 

49 

69 

89 

C0 2 (Q i ) a 

2 

3 

4 

j 



C0 2 (Q 2 ) a 

2 

4 

7 

( 92 

132 

170 

C0 2 (Q 3 ) a 

5 

10 

18 

' 




The numbering of normal coordinates is that used in the program — 
> > w 3 — and does not correspond to standard numbering. 
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APPENDIX 


SPECIAL PROBLEM - C0 2 


INTRODUCTION 


A subroutine and data set has been developed to generate the 
frequencies and intensities of Raman scattering lines of C0 2 in 
the 1300 cm- 1 region. This is to be incorporated into the general 
program for simulating the Raman spectrum of a combustion system 
as a function of temperature and number density, described 
earlier . 

The trace scattering contains the predominant features of the 
vibrational spectrum, since the depolarization ratio, p, is .027 
(ref. 2) for C0 2 , where 


P = 


31 . 

anis 

451, + 41 

tr anis 


I and I . are the trace and anisotropic scattering intensity, 
tr anis 

respectively. The trace cross-section for a transition from initial 
state, Ji, to final state, Jm, is given by 


x Ji , Jm 
da \ 
dft L 

/trace 


(2irw R ) 4 Yj _ 



a 



(A-l) 


where u. and u are the vibrational wave functions of the initial 
i m 

and final states, and the other quantities are as defined earlier. 
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The vibrational energy levels, and hence frequencies, of C0 2 
are not reproduced by the harmonic oscillator formulae given pre- 
viously, due to Fermi resonance between Vi and 2v 2 . it is 
first necessary therefore to determine the eigenvalues and eigen- 
functions of the perturbed Hamiltonian matrix. The unperturbed 
energy levels and the perturbation terms are calculated using param- 
eters from detailed analysis of infrared spectra (refs. 7, and 13 to 
18) . The resulting matrix is then diagonalized. In order to insure 
correct interpretation of the parameters, the energy eigenvalues, 
first obtained only for the nonrotating molecule, were compared 
to the calculated and the observed infrared values for 30 levels 
ranging from 667 to 7603 cm -1 . 

Two programs were used to perform the calculation: UNPERT, 

which generates the unperturbed frequencies and quantum numbers, 
and PERT, which reads the data from UNPERT, sets up and diagonalizes 
the energy matrix and outputs the energy eigenvalues. PERT was 
then - converted to'-ia LINES subroutine analogous to the rotational 
and diatomic vibrational LINES subroutines described before. 

This involved setting up the necessary data in terms of common 
blocks and arguments consistent with the basic program, solving 
the energy matrix for eigenvectors as well as eigenvalues, using 
the former to generate the cross-sections and the differences 
between the latter to generate the frequencies of the vibrational 
transitions. At the end of the subroutine the frequency and cross- 
section arrays are reordered according to increasing frequency 
in order to be consistent with the C parameter in the CNVLUT 
subroutine . 

Pure vibrational spectra were produced which did not include 
the dependence of the Fermi coupling terms on the rotational quantum 
number, J, nor a set of rotational lines for each vibrational 
transition, which forms the envelope of the vibrational band. A 
rough simulation of the band envelope was provided by giving a 
broad (3 cm’' 1 ) Lorentzian lineshape to each vibrational transition. 

The programs have subsequently been modified to include the 
rotational structure as follows: the unperturbed values of the 
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rotational constants, B^, were generated as a function of the 
vibrational quantum numbers, according to equation (9) , along 
with the E^. in program UNPERT. A loop and branch condition 
were added to the LINES subroutine to provide for explicit 
calculation of the rotational lines. This involved addition 
of the J dependence to the Fermi coupling terms and calcula- 
tion of the perturbed B v values by averaging over the eigen- 
vectors. These B v values are used to obtain the rotational 
contribution to the energy, which determines the rotational 
Boltzmann factor, and to obtain the rotational part of the 
frequency according to equations (7) and (8) . 

It is desirable to test the program at a rotational and vibra- 
tional temperature of 1565 °K, since an experimental spectrum at this 
temperature has been published (ref. 10). However, there are approx- 
imately 280 vibrational levels occupied at this temperature, each 
containing up to 145 rotational levels. Since it is not possible 
to load a two-dimensional array containing this many elements, the 
program was modified to provide for convolution over the line shape 
function after the lines for each set of vibrational levels for a 
given J value (i.e., rotational state) were generated. In this 
manner all the data is stored in the array of ordinate points 
(Y array as described earlier) , and the same frequency and 
intensity arrays are reused for each J value.- This is accom- 
plished by calling the CNVLUT subroutine from within the loop 
over J, which is the outer loop in the LINES subroutine. When 
j = nj - 1 (NJ is the maximum value of J) , control is transferred 
to the main program where the final call to subroutine CNVLUT is 
made. Convolution is then performed over the slit function, and 
the data is sent to subroutine SPKPLOT to produce a plotted spectrum. 

At present the unperturbed energy levels as calculated by an 
independent program are read in by this LINES subroutine. 
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THEORY 


Energy 

The expression for the unperturbed vibrational energies is 
essentially the same as that used before [see eq. (10) in main 
text] , except that the anharmonicity constants are the negative 
of those defined previously, and it is expressed more simply m 
terms of the anharmonic frequencies: 


G°(v 1 , v 2 , v 3 ) = Y 


k=l 


w k v k 


+ E E x., v. v + g 22 l : 

k=l j=l 3 K J K 


The w, are related to the harmonic frequencies wf by 
k K 


(A-2) 




= w. 


+ X 


kk kk 


I S 


x d . 
Dk j 


where the d, are the degeneracies of the k^* 1 normal coordinate, 

JC . 

i . e . , d 2 = 2, di = d 3 = 1. The term in in equation (A-2) was 

omitted in reference 1, and the degeneracy of Q 2 was not taken 
into account. 

The theory for the effect of Fermi resonance on the energies 
calculated by equation (a- 2) is essentially that of perturbation 
theory. Taylor, Benedict, and Strong (ref. 13) showed that the 
energy levels obtained from infrared spectra could not be fitted 
unless the Fermi coupling parameters were allowed to vary in a 
certain way with v 2 , and v 3 . Amat and Goldsmith (ref. 20) 

showed that the variation determined empirically by Taylor et al. 
is that given by perturbation theory when third-order terms in the 
transformed Hamiltonian are included. Amat pt al . - showed that the 
only nonvanishing matrix elements for this perturbation Hamiltonian 
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If v 2 


are the off-diagonal elements, (vj , v 2 , Z, v 3 | H 1 jv^ - 
+ 2, Z, v 3 y, which are also the only nonvanishing elements for 
the second-order case. Maes (ref. 21) determined the sum of these 
second- and third-order elements as a function of the quantum num- 
bers to be 


Vt(v 2 + 2) 2 - 
2 ' 


k 122 

— + \ l V X 

/2 


+ X 2 (v 2 + 2 ) 


+ X 3 (v 3 



(A- 3) 


Although coupling terms exist only between two unperturbed 
states at a time, any number of states can interact simultaneously. 
Symmetry requires that all interacting states have the same Z and 
v 3 values. The angular momentum quantum number, Z, must be 5 v 2 
and takes on alternate integer values, being even or odd in accord- 
ance with v 2 . 

The sets of interacting states, or "polyads," are shown system- 
atically for up to v 2 — 10 for v 3 = 0 in table 1. It is seen that 
for each value of v 2 there is one polyad each whose degree n is 
given by 

v 2 - £ 0 

n = 1 , 2 , 3 , . . . 9 + 1 


where Zq minimum value of Z: 


Z 0 = 0 for v 2 even; 

£o = l f° r v 2 odd. 

Referring to the right side of table Al, it can also be seen that 
for each degree n there is one polyad of that degree for every Z 
value; Z - 0 , 1, 2 ,...,& , where 

IIlclX 
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Table Al 


Organization of polyads in terms of 
unperturbed states (v lf V2 , V3) . 

























































2 (n - 1) . 


„ , , max 

& (n) = v 2 

max z - 


Within a given polyad, if we consider the levels as given in 
the table from left to right, i.e., starting with v 1 - 0 and 
increasing it consecutively by one and decreasing v 2 corres- 
pondingly by two, the energy matrix takes on " tridiagonal " form: 


E* W 2 0 0 0 

W 2 E° W 3 0 0 

[<u*|H|up] = 0 W 3 E» W 4 0 


0 0 W 4 E 4 


( A-4 ) 


The v x and v 2 values may be set up in this order in terms of 
l r n, and a running index, k : 


V! = k - 1 

v 2 = 2 (n - k) + & 


k = 0, 


1 , 


2 , 


n 


(A-5) 


This index is shown in table Al. The E£ in equation (A-4) are thus 
the unperturbed energies calculated from equation (A-2) , and the W k 
are the coupling terms calculated from equation (A- 3) , with the k 
subscripts determining the V]_ and v 2 quantum numbers. Diagon- 
alization of this matrix which yields the correct energy eigenvalues 
for the polyad characterized by v 3 , &, and n, is carried out by 

the QL algorithm (ref. 22) , using the subroutine TQL2 , which is 
available in the CDC math library. 


Frequency and Cross-section 

The selection rules for the Raman scattering transitions are 
determined by the symmetry of the polarizability operator and the 
fact that we are interested in the 1300 cm” 1 region which corresponds 
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to transitions, Av : = 1 or Av 2 = 2, for the unperturbed levels. 
With the Fermi resonance the Vj and v 2 quantum numbers are no 
longer valid, and each state in the polyad becomes a mixture of 
all the unperturbed states in the polyad. Thus, if one state in 
a polyad can undergo a transition to a state in a second polyad, 
all the states in the first polyad can undergo a transition to 
each of the states in the second. Since Qi and Q| are both 
totally symmetric, symmetry requires that A& — AV 3 = 0 for the 
transitions. It can be seen from table Al that all the transitions 
for which Av! = 1, Av 2 = 0, and Av 2 = 2, Av x = 0 for the unperturbed 
states are contained in the conditions An = 1, A£ = 0. Hence the 
selection rules for the perturbed states are: 

A£ - A v 3 = 0 

and 

An = 1 . 


The transition frequencies are determined simply by subtracting 
the energy of the initial state from that of the final state. The 
expression for the transition cross-section is obtained from equa- 
tion (A- 1) , using 


u. 

1 




(A- 6 ) 


where the u°. are the unperturbed wave functions and c.^, 
j = is the i eigenvector of the matrix m equation 
(A- 4 ) . The relevant terms in the expansion of the polariz- 
ability operator a are 


a 


,8 a 

9Qi 


Qi + 


1 ' 2 

o 8 Q 2 


9 aX qI 


) 
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where a represents the trace of the polarizability tensor as 
defined before, equation (19) in main text. Representing the 
by (vi , yj? ), since it is understood that the H and v 3 quantum 
numbers of the initial and final states are equal, we have 


n 


/\ 

/u a u.'y 
v m 


E 

k=l 


n-1 

E 

j=i 


'km ji 




k \ 

v 2 ) 


a 


(v? / 


v 2 )> 


n n-1' 


= E E 

k=l j=l 


c. c . . «k - 1, 2 [n - k] + A) 
km ji N 


(A-7) 


x 



(j - 1, 2 [n - 1 


j] + O 


where we now express the v 3 and V 2 quantum numbers in terms of 
the running index defined in equation (A- 5) , bearing in mind that the 
degree of the polyad to which the initial state belongs is n — 1. 

The sums in equation (A-7) may be separated into two types of terms 
according to the integrals over the normal coordinates in the unper- 
turbed wave functions. We have (refs. 24, 25) 


_ i i \ 3a \ 

v i + i|qKi> = 

is the integral evaluated for v x = 0, and 


3a\ 

9Qij 


< 


where (a{)' 


V 


h(v x + 1) 
8tt 2 v i 


= <ai>° /v~i + T (A-8) 


3 2 a 
9Q 2 


\ <v 2 + 2 | Q 2 | v 2 > - - 

) 0 9Q ^o 


h 


16tt 2 v 2 


V(v 2 + 2) 2 - £ 2 


(A-9) 


<a 2 > 


V (v 2 + 2) 2 - z 2 
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where <(a 2 >° is the integral evaluated for v 2 = £ = 0. It can be 
seen that the first type of term occurs for j = k - 1 in equation 
(A- 7) and the second type for j = k. Accordingly, equation (A-7) 
reduces to 


n 


<jj | a | u .y = 
x m 1 r 


k=2 


c km c k-l,i 


- 


Ql |k - 2> 


n-1 


£ c k- °M <2|n ■ kl * l| !s) 


2 ™ X Q 2 | 2(n - 1 - k) + 


where the u° 1 s are now represented by only one quantum number — 
that which is different for the two states — v x for the terms in 
the first sum,- v 2 for those in the second. The sums over k are 

altered due to the fact that the sum over j goes from 1 to n - 1; 

thus, there is no j = k term for k = n, and no j =-k - 1 term for 
k = 1. Replacing the dummy index k by k + 1 in the first term, 
we see that it is given by equation (A-8) with vj = k - 1. The 

second term is given by equation (A-9) with v 2 + 2 = 2(n - k) + &. 

Thus we have 


n-1 


< u J“i u i> = <“i>° k ^ i c ki 


c, ,, /k 
k+1 ,m 


+ c 


km 


<«2> U 

<a,)° 


V (n - k) (n - k + i) 


(A- 10) 


This quantity is then squared and multiplied by the appropriate 
factors according to equation (A-l) in order to give the cross- 
section for the transition. 
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Rotational Structure 


Maes has shown that the Fermi coupling term has significant 
J dependence for rotational levels with J > 30 (ref. 21) . These 
levels will have large populations at the temperatures considered 
in combustion systems. He has determined the functional form of 
this dependence to be 


Vt(v 2 + 2) 2 - & 2 3 v ! 

<$J(J + 1) 5 


( A- 1 1 ) 


This term is accordingly added to the W^. in equation (A- 4) , as 
given by equation^ (A- 3) , when the rotational structure is con- 
sidered. The eigenvalues’ and eigenvectors of equation (A-4) will 
then be different for each J value. 

The rotational energies are given by equation (8) in the main 
text, except that now the B v values msut be averaged over the wave 
function for the perturbed states : 

b ;= e (4>* (a-12 

k=l 


where the J dependence of the eigenvectors and, hence, of the 
perturbed B v values is now reflected. 

The rotational selection rules for the trace scattering are 
given by A J = 0 ; hence, f the rotational frequencies are given by 

J J J ' 

coj = (B m - B ± ) J(J + 1) (A_ 

where i and m - represent the vibrational quantum numbers of the 
initial and final states as before. 

The J dependence of the rotation-vibration cross-sections 
derives from the J dependence of the rotational energies and 
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hence, of the rotational Boltzmann factor, and the J dependence 
of the eigenvectors as reflected in the transition polarizability 
given by equation (A-10) . There is no rotational line strength, 
S rt , for the trace scattering. 

DATA AND RESULTS 

The parameters used in the calculations are those determined 

by Suzuki (ref. 18) (set v) and are given in table A2 . They were 

obtained by a simultaneous fit of the parameters indicated to 33 

vibrational band origins and 31 rotational constant differences , 

(B - B ) , observed from infrared spectra. The ground state 
v o 

rotational constant f B^ = .39021 cm, was that used in reference 1, 
first obtained by Courtoy (ref. 7) and verified by a more accurate 
measurement (ref. 17) . 

The comparison of those energy levels obtained by the programs 
described above, those observed from infrared spectra, and those ca 
culated by Suzuki is shown in table A3. It can be;, seen that the 
agreement is reasonable. Further, it can be seen that the agree- 
ment is remarkably good between the calculated energies and those 
observed energies to which the parameters were not fitted (those 
for which there are no entries in column 2) . 

In addition to the parameters listed in table A3, the polar- 
izability ratio, (®2y°/*\ a 1 /*° r an< 3 the value, <^i)> 0 , in equation 
(A-10) are needed in order to calculate; the spectrum. The value of 
the polarizability ratio was determined from the measurement of 
the ratio of the Raman intensity of the two fundamental bands 
carried out by Howard- Lock and Stoicheff (ref. 25) . The ratio 
was obtained using equation (1) of that reference, but using the 
energy parameters of Suzuki rather than those of Courtoy as did 
Howard- Lock et al. 

The value of (ax)° was determined from the absolute cross- 
section measurement of Penney et al. (ref. 12) using the equation 













Table A3. Calculated and observed energy levels 
for COo, cm” 1 (concluded). 


Observed 

Calculated 
(Suzuki, ref. 18) 

Calculated 
(present work) 

4977. 81 d 

4977.92 

4977.75 

5061. 78 a 

5062.07 

5062.21 

5099. 61 d 

5099.64 

5099.69 

5316. 09 d 

5316.18 

5316.21 

5960. 08 d 

5959.42 

5959.49 

6075.93 d 

6075.52 

6075.69 

6227. 88 d 

6227.76 

6227.86 

6347. 81 d 

6348.25 

6348.35 

6503. 05 d 

6503.12 

6503.30 

6688. 54 d 

6687.57 

6687.59 

6863. 91 d 

6863.77 

6863.94 

6972. 49 d 

6972.70 

6972.63 

7024.03 d 

7024 . 50 

7024.68 

7204. 22 d 

7204.01 

7204.23 

7602. 85 d 

7603.12 

7603.11 


Reference 17. 

Reference 14. 

Reference 5 as quoted in reference 7 
Reference 7 . 








zz 


3 , ~ \ 4 V-0 

3 - 4p ( 27 TW r ) n T=298 


«<*!>> 


1 + 


<ot 2 > C 

< a i> C 


where p is the depolarization rato and n T=298 t ^ ie re l at: *- ve 

population of the ground state at room temperature. All values 
used in the equation were those of Penney et al. except for the 
polarizability ratio, which was determined as described above. 

It was assumed that the absolute cross-section measurement was 
over a frequency region including both the fundamental bands. 

The values determined were 

<^ 2 >° 

= -.126663 , 

<ai>° 


<a !>° = 1.2215 x 10 26 cm 3 


The spectra with and without rotational structure are shown in 
figures Al and A2 , - respectively . The instrument parameters and 
number density, which determine only the scale of the spectrum 
and not its relative shape, are the same as those used in the 
main text. The remaining! input data are given in table A4 . It 
can be seen that explicit calculation of the rotational structure 
makes a substantial difference in the spectrum. 

Figure A3 shows the experimental C0 2 spectrum taken by Lapp at 
1565°K (ref. 19) . It can be seen that the calculated vibration- 
rotation spectrum is in good agreement with this spectrum. 
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Figure A2 . Vibrat 






















(03 3 0)-(l3 3 0)+ (05 3 0) (00°0)-(l0°0)+(02°0) 



Figure A3. Experimental C0 2 Stokes vibrational Raman spectrum for C0 2 seeded into a 
stoichiometric H-air flame at about 1565°K, and at a partial pressure of 
about 1/3 atm. The data were obtained through use of a l.lw, 488 nm 
argon ion laser source and a 3/4-meter double monochromator with a tri- 
angular spectral slit function full width-half maximum value of 0.163 nm. 
The notation for the indicated transitions, (vj , v 2 , £, V3), corresponds 
to the three fundamental vibrational quantum numbers and to the vibra- 
tional angular momentum £ of the v 2 -bending mode. 
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